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From first-principles approaches, we illustrate that the current-induced forces and the selection 
rule for inelastic effects are highly relevant to the current density in an asymmetric molecular 
junction. The curved flow of current streamline around the asymmetric molecule may induce a net 
torque, which tends to rotate the benzene molecule, similar to the way a stream of water rotates 
a waterwheel. Thus, the Pt/benzene junction offers a practical system in the exploration of the 
possibility of atomic-scale motors. We also enumerate examples to show that the use of selection 
rule can lead to misjudgement of the importance of normal modes in the inelastic profiles when the 
detailed information about the current density is not considered. 



Electron transport in a single-molecule molecular junc- 
tion, where the molecule is sandwiched between two elec- 
trodes, has been investigated extensively in the pursuit of 
extreme device miniaturization [H-Q. A major concern 
for single-molecule junctions is the fundamental proper- 
ties of current-induced effects due to nonequilibrium elec- 
tron transport at finite biases. These effects are efficient 
tools to explore the single-molecule signatures from quan- 
tum mechanical perspectives, and thus are important and 
interesting from theoretical and experimental points of 
view. For example, the question on whether the current- 
induced forces on atoms are conserved has been raised @ . 
Through molecular dynamics simulations, it has been re- 
cently reported that the current-induced forces are gen- 
erally not conserved @, 0|- The authors show that it 
is possible to construct atomic-scale systems where the 
current-induced forces can be used to rotate the atoms. 
Therefore, the current-induced forces are not conserved 
due to the nonzero net work done by the current-induced 
forces 0- 

The current-induced effects with nuclear degrees of 
freedom of molecules can enrich the understanding of 
the interactions between electrons and molecular vibra- 
tions. In the hydrogen junctions, the mode identification 
and the selection rule for modes that significantly con- 
tribute to IETS are related to the large component of 
vibration along the direction of electron transport (z- 
direction) Due to the short length of the junction, 
electrons are transported quasi-ballistically in the molec- 
ular junction 0. Consequently, both the interaction be- 
tween electrons and local ions (electron-vibration inter- 
action) and the induced local heating need to be consid- 
ered. 0-tii- Electron- vibration interaction reveals the 
junction signatures through the inelastic electron tunnel- 
ing spectroscopy (IETS) [3, EH, while local temperature 
resulting from local heating can be obtained by equili- 
brating the heat generation in the junction with heat 
dissipation to the electrodes. Recent report shows that 
local heating leads to a stronger suppression of the cur- 
rent at longitudinal vibrational modes for a 4-A1 atomic 
wire flij ]. 

Moreover, Kiguchi et al. have measured a high con- 



ductance in a single-molecule junction where the benzene 
molecule is connected directly to the platinum electrodes 
(Pt/benzene junction) [l7|- This demonstrates the feasi- 
bility of a carbon-metal link. The authors have concluded 
that the benzene molecule forms a direct bond with the 
electrodes. This is based on the analysis of the conduc- 
tance histogram and IETS, where the vibrational modes 
in the inelastic profiles have been verified by isotope sub- 
stitution. The relaxed configuration of the Pt/benzene 
junction, as we shall show later, loses mirror symmetry. 
The highly tilted benzene molecule causes the stream- 
line flow of the current to curve considerably to one side 
of the benzene ring. This could cause the unbalanced 
current-induced forces, which tend to rotate the benzene 
molecule. Moreover, the mode selection rule for IETS 
in the highly tilted benzene molecule junction becomes 
intricate. The selection rule based on the component of 
vibrations along the z-direction turns out to be inappro- 
priate. Determining the importance the contribution of 
normal modes to IETS based only on the z-component of 
vibrations will lead to misjudgement. This suggests the 
necessity of looking into the detailed current density, es- 
pecially in highly asymmetric molecular junctions, such 
as the Pt/benzene junction. 

In this Letter, we demonstrate that the current- 
induced effects are highly relevant to the details of the 
current density. This is particularly important in the 
Pt/benzene junction because of its highly asymmetric 
configuration. From first-principles approaches, we in- 
vestigate the impact of asymmetric current streamline 
flow on the current-induced forces and the selection rule 
for IETS in the Pt/benzene junction. It should be noted 
that our theoretical method has an advantage over the 
DFT+nonequilibrium Green function method since the 
latter does not offer information about current density. 
We show that the curved current streamline flow around 
the molecule may induce a net torque, which tends to 
rotate the benzene molecule. To compare with the ex- 
perimental IETS data, we calculate the inelastic profile 
of conductance (dl/dV) and the derivative of conduc- 
tance (d 2 I/dV 2 ) with and without local heating. Coun- 
terexamples are provided to show that the selection rule 
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for important normal modes based on the component of 
vibrations along the z-direction is inadequate. 

We start with a brief introduction to the theories which 
lead to the results of our calculations. We study the 
current-induced forces and inelastic profiles in the frame- 
work of density functional theory (DFT) in scattering ap- 
proaches [HI • The wavefunctions $^ (r) are calculated 
by solving the Lippmann-Schwingcr equation iteratively 
until self-consistency is obtained. 

The force F, acting on a given atom at position R 
exerted by the nonequilibrium current, is given by the 
Hellmann-Feynman type of theorem which has been de- 
veloped in Ref. [H. 
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where the first term on the right hand side of Eq. [T] is 
the Hellmann-Feynman contribution to the force due to 
localized electronic states The second term is the 

contribution to the force from the continuum states |"0a}- 
The inelastic current and local heating are calculated 
using the first-order perturbation theory based on the 
second-quantized formalism. By applying the field oper- 
ator with the wavefunctions obtained in DFT, the many- 
body Hamiltonian of the system is H = H e i + H v n, + 
Hei-vib-, where H e i is the electronic part under adiabatic 
approximations; H v n, is the ionic part of the Hamilto- 
nian, which can be casted into a set of independent simple 
harmonic oscillators via normal coordinates; and H e i- V n, 
represents the electron- vibration interactions and has the 
form 
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where a and /3 refers to either the left (L) or the right 
(R) electrodes; Mi is the mass of the i-th atom; A^j is a 
canonical transformation between normal and cartesian 
coordinates satisfying Ai^jAi^j' = 8j,j>] uij are the 

normal mode frequencies; bj is the annihilation operator 
for phonons corresponding to the j-th normal mode; and 
the annihilation operator for electrons. The cou- 
pling constant Jj^'^f between electrons and the vibration 
of the i-th atom in /i-th (/x = x, y, z) component, i.e., 

J&% =J dr J rfKiiI^K^n^^R^iK/r)] 

(3) 

where V ps (r, Hi) is the pseudopotential representing the 
interaction between the electron at r and the i-th ion at 

Ri- 

The current in the presence of electron- vibration scat- 
tering within the first-order correction is, flij ]. 

I = ^ / dE[(f§ fk) {B R B L )]t(E), (4) 



where = l/[cxp((£ - VL{R))/(k B T L{R) ))+l] is 

the Fermi-Dirac distribution function for the left (right) 
electrode, ks is the Boltzmann constant, fJ-L(R) the 
chemical potential in left(right) electrode, TWm is the 
temperature in the left(right) electrode, and t(E) = 
z€jdRj dK||(*f*V*f - V#§**|) is the transmis- 
sion function of electrons with the energy E from the 
right electrode. In the absence of electron-vibration in- 
teractions, Eq. (Ql returns to the Landauer-Biittiker cur- 
rent formula. The normalized parameters B a (a = L, R) 
from the inelastic scattering processes are obtained by 
the following equations: 
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where {a, /?} = {L, R} and a ^ (3. The parameter 
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where 5 = 0(1) represents the process of phonon absorp- 
tion(emission) , and rij V = l/[cyLj>(hkjj V /kBT w ) — 1] de- 
notes the ensemble averages of phonon states where T w 
is the local temperature. (Bj" k ) refers to the ensemble 
average over phonon states. 

To incorporate local heating in our calculations, we 
compute the total thermal power generated in the junc- 
tion P via electron- vibration interactions using flij , 
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where W"^' 2 ^ 1 ' is the power calculated up to the first 
order from the Fermi golden rule, corresponding to re- 
laxation (excitation) of the vibrational modes [10( . 

We consider successively the effects of the curved flow 
of current streamline on the current-induced forces and 
inelastic profiles in a highly asymmetric single-molecule 
junction. First, we compare the current density at 0.1 V 
with the forces on each atom in the Pt/bcnzcne junction, 
where the current-induced forces are calculated by sub- 
tracting the total forces at 0.1 V and V, as shown in 
Fig. [IJa). We observe that the benzene is subjected to 
comprcssional forces in both the horizontal and vertical 
directions. This feature is consistent with the current- 
induced forces for the conventional Au/benzene junction, 
where the benzene molecule lies on the same plane with 
both electrodes and exhibits mirror symmetry [20| . 

In the Pt/benzene junction, the horizontal compres- 
sional force on the two outermost platinum atoms is sig- 
nificantly stronger than the vertical force because the 
current flows mostly through these two atoms. Due to 
the tilted benzene ring with respect to the y — z plane, 
the net force is considerably stronger in the negative y 
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mode has also been observed in the IETS but has not 



FIG. 1: (Color online), (a) 3D force vector plot (blue line) 
of the benzene molecule (red) connected to platinum atoms 
(purple) with Vb = 0.1 V on top of the current density vector 
plot (black line). The magnitude of the current density varies 
in orders of magnitude. The graph is plotted in log scale 
to help the visualization of the smaller current density. The 
separation of Pt-electrodes (jellium modeled, r s ~ 3) is fixed 
at 9.626 a.u. (b) Projection of the force vectors on the plane 
of the benzene molecule. The projected net torque around the 
center ofbenzene ring is (0.0013, -1.637 x 1CT 7 , -0.0016) nN- 



and negative z directions. The curved current streamline 
flow around the molecule produces a net torque around 
the center of the benzene ring. The net torque [around 
(0.010, -0.077, 0.005) nN-A] is stronger in the negative 
2/-dircction. This means that the force induced by the 
curved current stream line flow could lead to a clockwise 
rotation in the x — z plane. We also project the force 
vectors onto the plane of the benzene ring, as shown 
in Fig. [Tf b) . We observe that the benzene ring expe- 
riences a net torque which tends to rotate the molecule 
clockwisely. Thus, the highly asymmetric single-molecule 
junction, such as the Pt/benzene junction, offers a prac- 
tical system in the exploration of the possibility that cur- 
rent could be used to drive atomic-scale motors. 

Second, we show that the current streamline flow is 
crucial to explain the selection rule for normal modes 
shown in the IETS. Without looking into the details of 
current density, the selection only based on the current 
could lead to misjudgement with regard to the impor- 
tance of normal modes for IETS in the Pt/benzene junc- 
tion and other highly asymmetric molecular junctions. 

In Fig.[2ja), four modes denoted as (I), (II), (III) and 
(IV) out of 42 vibrational modes with normal mode ener- 
gies around 42 eV, 64 eV, 80 eV, and 93 eV, respectively, 
are shown. Mode (I) is acoustic-like and contributes to 
IETS at Vb = 42 mV, which has been experimentally 
observed (l7j . This mode corresponds to a vibration of 
the benzene molecule as a whole, with a large compo- 
nent of motion along the current. In other words, the 
center of mass of the benzene molecule moves along the 
^-direction. Mode (II) is a longitudinal mode where the 
middle part vibrates in the direction opposite to those 
of the top and the bottom parts. Mode (III) is a trans- 
verse mode where the left side of the carbon ring vibrates 
downwards while the right-side vibrates upwards. This 



been discussed in Ref. [17[. Mode (IV) is another longi- 
tudinal mode. However, only the second-nearest neigh- 
bor vibrates in the same direction. The center of mass of 
the benzene molecule remains fixed. 

To include the effects of local heating, we obtain 
the rate of thermal energy dissipated to the bulk elec- 
trodes via phonon-phonon interactions using the weak- 



link model [21|. The stiffness of the benzene molecule is 
estimated using the total energy calculation, which gives 
K = 6.67302 x 10~ 4 cV /a§, where a is the Bohr ra- 
dius. The effective local temperature T w is achieved when 
heat generation in the nanostructure and heat dissipa- 
tion into the bulk electrodes are balanced • Fig. ^h) 
shows the effective local temperature T w vs. Vb in a 
Pt/benzene junction, where the temperatures of bulk Pt 
electrodes are set to K. We observe that the local tem- 
perature increases rapidly when the bias is larger than 
Vb= 42 mV, where the electrons have enough energy to 
excite mode (I), the lowest energy mode vibrating along 
the z-direction. 

Subsequently, we calculate the inelastic profile of con- 
ductance (dl/dV) [upper panel in Fig. [^c)] and the ab- 
solute value of the differential conductance (d 2 I/dV 2 ) 
[lower panel in Fig. [2jc)], with and without the presence 
of local heating, as a function of bias voltage. The zero- 
bias conductance (approximately 0.58 Go) is in agree- 
ment with the large conductance observed in the ex- 
periments 17J. This is in sharp contrast to the small 



conductance which has been found in the conventional 
junction, where the benzene molecule is connected to 
the gold electrodes via the sulfur atoms. When local 
heating is included, we observe that the inelastic fea- 
tures are enhanced due to increased number of local 
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FIG. 2: (Color online), (a) Four vibrational modes corre- 
sponding to energies (1)42 eV, (II) 64 eV, (III) 80 eV, and 
(IV) 93 eV. (b) Local temperature T w as a function of the ap- 
plied bias Vb at zero electrode temperature, (c) Conductance 
and the derivative of conductance as a function of the applied 
bias Vb with local heating ( red and dashed) and without lo- 
cal heating(black and solid). Temperatures of electrodes are 
set to zero. 
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FIG. 3: (Color online), (a) Schematic of Mode (III) and 
(b) Mode (II) with vibrational direction (blue line) on top 
of the current density vector plot (black line) of the ben- 
zene molecule (red) connected to the platinum atoms (purple) 
sandwiched between two Pt electrodes with a bias voltage of 
V B = 0.1 V. 

phonons. Figure |5Jc) shows that large jumps in IETS 
correspond to three normal modes [Modes (I), (III), and 
(IV)]. Modes (I) and (IV) clearly show a large compo- 
nent of vibration along the direction of electron transport 
(i.e., roughly along the line connecting two Pt atoms), as 
shown in Fig. [Ha). The large contributions of these two 
modes can be explained by the selection rule based on 
the current: the strength of inelastic effects are related 
to the component of vibration along the z-direction Q . 
However, Mode (III) is an exception to the selection rule 
based on the current. Fig.^a) indicates that Mode (III) 
is a transverse mode. The motion is mostly perpendic- 
ular to the line connecting the two Pt atoms, and thus, 
its contribution to the inelastic effects is expected to be 
minimal. The large contribution of Mode (III) to the 
IETS has posed a perplexing problem regarding the ap- 
plicability of the selection rule based on the current. 

To settle this problem, we visualize the current density 
and compare the electron current streamline flow with 
the motion of normal modes. It should be noted that the 
center of the mass (CM) of the benzene molecule does not 
lie in the line connecting the two platinum atoms. Under 
such highly asymmetric configuration, the current tends 
to flow through the bottom half of the benzene molecule, 
as shown in Fig. [3J For Mode (III), the carbon ring near 
the left electrode vibrates upward while the carbon-ring 
near the right electrode vibrates downward. Fig. [3ta) 
clearly shows that the carbon-ring vibrates along the 
curved electron current streamline flow. This justifies 
the importance of Mode (III) to the IETS. We illustrate 
another example which shows the importance of the cur- 
rent density on the selection rule for inelastic effects. Fig- 
ure [He) displays that Mode (II) has a large component 
of motion along the z-direction, but the corresponding 
inelastic profile does not occur, as shown in Fig. ^c). 
Figurc^b) shows that the motion of Mode (II) is mostly 
transverse to the curved current streamline flow. Thus, 
this mode is unimportant to the IETS. 

In conclusion, we have shown that the current-induced 
forces and inelastic profiles are highly relevant to the de- 
tails of current density, especially in highly asymmetric 
molecular junctions. The asymmetric current streamline 



flows may cause net torque, which tends to rotate the 
molecule, in a manner similar to a stream of water ro- 
tates a waterwheel. Thus, the highly asymmetric single- 
molecule junctions offer practical systems to explore the 
possibility of an atomic-scale motor driven by the non- 
equilibrium electron transport. Current density is also 
important in the explanation of the mode selection rule 
for inelastic effects. In the absence of detailed informa- 
tion on the current density, we enumerate examples to 
show that the importance of normal modes on the in- 
clastic effects could be misjudged. 
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